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We investigate a particular two-point function of the n-copy Ising model. That is, the corre- 
lation function (e(r)T(O)) involving the energy field and the branch-point twist field. The latter 
is associated to the symmetry of the theory under cyclic permutations of its copies. We use a 
form factor expansion to obtain an exact integral representation of (e(r)T(O)) and find its com- 
plete short distance expansion. This allows us to identify all the fields contributing in the short 
distance massive OPE of the correlation function under examination, and fix their expectation 
values, conformal structure constants and massive corrections thereof. Most contributions are 
given by the composite field : eT ■ and its derivatives. We find all non-vanishing form factors of 
this latter operator. 
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1 Introduction 



The idea of associating a local field to a branch point in a CFT on a Riemann surface was first 
introduced in [I], and deepened in [2]. It became then natural in the context of the evaluation of 
bi-partite entanglement entropy for two-dimensional integrable QFTs (see i.e. [3] for a review), 
and for non-integrable ones These models are usually quantized on a line, and then evolved 
in euclidean time r = it. The idea is to divide the system at r = into two regions, and then 
consider n copies of the original model. These "replicas" are connected cyclically along a branch 
cut which covers one of the two regions (the choice is indeed arbitrary). The boundary points 
of the cut are branch points of degree n, and a field T is associated to them, related to the Z n 
symmetry of the global model. These fields are called twist fields, and are local with respect 
to the hamiltonian density of the n-copy model h n (x) = Y^j=i^j( x )i where hj(x) is evaluated 
in the j-th. sheet of the n-root map which connects the Riemann surface with M 2 . Consider for 
simplicity Aio the manifold with a single branch point in x = and the disconnected n 
replicas, then the action of the twist field on the expectation values of the theory is defined as 

{■■■)Mo = pP[ • l 1 ) 

The counterpart of the twist field in the ultraviolet conformal theory is a spinless operator of 
conformal dimension [2] 

c / 1' 



A T = A T =-\n--), (2) 

In a recent paper [5] involving the author, it was first observed that a function that generalizes 
A 7- to massive theories, known in literature as A-function [6] has all the right properties to 
be a c-function of the theory. That is it has the same qualitative features of Zamolodchikov's 
c- function [7], even if it is different from it. This observation was put on a more solid basis in 
[8], where the correlation function (0(i)T(O)) between the perturbing field and the twist field 
in general 1 + 1-dimensional integrable theories was considered. In particular it was noticed 
that it can be obtained exactly for the n-copy thermally perturbed Ising model, where <j> = e, 
that is the energy field of conformal dimension A = A = 1/2. This gives us the opportunity 
of investigating the massive OPE for the operator product e(r)T(0) in great detail. This is the 
main issue we want to address in this paper. 

The present work is organized as follows: after Section [21 that is devoted to describe the model 
and set the notation, the analysis of the two point function (e(r)T(O)) and the identification of 
the operator : eT ■ is carried out in Section [3l The two-particle form factor of this composite 
operator is found in Section [U and higher-particle form factors in Section [5j In Section [6] we 
collect all our main conclusions, and discuss open problems. 



2 The Ising model and twist field form factors 

The two dimensional thermally perturbed Ising model can be described by a relativistic quantum 
field theory on a line with action 

S E = i f d 2 x (tP(d x - d t )i> - iP(d x + d t )i> - mjn/>) , (3) 

which can be obtained by perturbing the critical Ising model with the energy operator and 
introducing a mass scale m which is related to the temperature. While in the conformal theory 
this field can be taken to be the product of two fermionic fields, in the off critical model, due to 
the presence of m, there is an additional mixing with the identity operator, so that 

e = a^tp + bml with a,6eR\{0}. (4) 
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In ((3|) i/j and ip are the two real components of a Majorana spinor, and they can be represented 
as 

i/}(t,x) = [ dee^a(9)e m{ixsinhe - TCOshe) + a^9)e- m{ixsinhe - TCOshe) } 

\l Air I 

(5) 

" N m(ix sinh 6— r cosh 9) f / q\ — m(ix sinh 9— r cosh 9) - 



if>(T,x) = -i\ — / d9e-2{a{0)e m{ - lxsmb °- TCOStl0 > - a T( ) e -^smii y -rcos ±iyj | ) 

V 47T J 

where 6 is the rapidity, and the mode operators a(9) and a^(9) satisfy the canonical commutation 
relations 

{a(e),a f (e')} = 6(9-0') {a(9), a{9')} = {a^{9), a^{9')} = 0. (6) 

The Hilbert space of this model is the Fock space built over the postulated vacuum |0) and the 
algebra © . The action of the operators a and a) is 

a(0)|O) = (O|at(0) = O at(0)|O> = |0> (0| a{9) = (0\ , (7) 

and a general vector is defined by multiple applications of ([7|) 

\9 1 ,9 2 ,...,9 n ) =a\9 l ) ( J(9 2 )...a\9 n )\U) 1 (8) 

where a basis is chosen setting 9\ > 9i > ... > 9 n . 

The fields of the theory act as linear operators on the Hilbert space. For an operator O the 
general matrix element 

if (0 1; ...A) = <O|0(O)|0i,...,0 fc > in , (9) 

is known in literature as a form factoiE- These quantities are particularly useful in the context 
of integrable theories, and they provide a non-perturbative expansion for two-point functions. 
Indeed by the introduction of the resolution of the identity over the basis (|8|) it is possible to 
express the correlator between two operators 0\ and 2 aqll 

(O x (r)0 2 m = g / ^^if ^i, • • • , k ) (F^{9 X , 9 k ))* e — Ei=x«** 

k=1 9 1 >9 2 >->9 k 

(10) 

The convergence of this series is remarkably fast for most integrable models, and leads to a good 
approximation just considering the first few terms. The aforementioned constants a and 6 in (J3]) 
are fixed by the normalization 

e(x) = 2vr (ml+ : tp(x)tp(x) :) , (11) 

where aa^ =: aa^ : +1 is the definition of normal ordering in this case. The two-particle form 
factor can be easily extracted using (JSJ) 

-i 
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F^9 U 9 2 ) = (O|e(O)|0i0 2 ) = -i? / d^eV (0| a( ?? )a( ( /))a t (ei)at(^) |0) 



. t/i — v 2 
-im sinh ■ 



(12) 



2 ' 

where the Wick theorem for fermion algebrae is used to define contractions. The normalization 
(|lip is chosen in light of the relation between e and the trace of the stress energy tensor Q(x) = 
2Trme(x) [9], to give F 2 e (z7r) = 27rm 2 . 



*in general form factors depend on a complete set of quantum numbers. The one-copy Ising model though 
has only one kind of "particle" in its spectrum, so that a form factor needs no other specifications than rapidities 
of incoming particles, whereas for the n-copy theory there is an additional number that has to be specified that 
labels the copy the "particles" belong to. 

^throughout this paper the following notation will be adopted 0(r = 0, x) = 0(x) 
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2.1 Form factors construction for T 

The twist field is a field in the general meaning of the term: is an object for which expectation 
values and exchange relations with all the other fields in the theory can be defined. In particular 
its action on the fermion fields is 

<Pi (t, y)T(r, x) = T{t, x)ip i+1 (r, y) x 1 > y 1 , 

^i{T,y)T(T,x) = T(T,x)tpi(T,y) x 1 <y 1 , (13) 

where ipi is the fermionic field of the i th copy, and the same relations hold for ^j. 
These non-trivial relations modify Watson's equations 

^ r|ni - nin ' I+1 -(- .-A, 0*1, • • • ) = sS m (^ i+ i)^ r| - remrei -(- • • , e*u Oi,...), (14) 
Fp nina - n "(e 1 + 2m t ...,e k ) = FP n2 - m+ \e 2 ,...,e k ,e 1 ), (is) 

where particles' indexes rii label the copy they belong to, 9a+i = B% — 9i+i and, for the Ising 
model, S^l +1 (6 ii+1 ) = 

The pole structure is also modified by (|13|) which introduces an extra kinematic equation as 
follows 

Res F[l n 2 ^ ni - nk (e Q + i7T,e ,e 1 ,...,e k ) = iF^ ni - nk (9 h ...,9 k ). (16) 

9q = 9q 

Res F™ ni+1)nt " n \9o + ™,9 ,0u---,0k) = -i {[ S^ rij (9 0j )F k Tl " k (9 1 , . . . , 9 k ^7) 
9q = 9o j=i 

These equations connect form factors of different orders, and can be used to check the validity 
of any form factors associated to the twist field. 

The two particle form factor for the Ising model was found in [10] . and can be expressed as 

F T|iW_ «(T)cos(^) sinh(^) 
2 U " sinh(^) sinh(^)' ^ 

and the exact expression of (T) was also given in that paper. Moreover higher particle form 
factors were also extracted, and they showed a Pfaffian structure. As T is an even operator 
under the Z2 symmetry of the Ising model only even particle form factors are non-vanishing, 
and they have the form 

F^ n - 1 (9 l ,...,9 2k ) = (T)Pi(K), (19) 
where Pf(if ) 2 = det(K) is a Pfaffian, and K is an anti-symmetric 2k x 2k matrix, with entries 

K ij ~ JY) ' *• ' 

3 An important OPE and composite twist fields 

In [8] the authors were interested in the computation of the delta-function of the twist field 

A(r,= -2i. ^hfk — (e H- <21> 
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where O is the trace of the energy-momentum tensor. There, evidence was given that this 
function has all the right properties to be a c-function. In light of the relation between and 
e mentioned in Section [2] we will henceforth refer directly to this second operator. Indeed for 
e only the two-particle form factor is non-vanishing and it is possible to find an exact integral 
representation for the correlation function^ 



oo 



, , Sr ^. NV , w _ v (T)m tt f K (2mrcosh~)sinh^-sinh§ t s 

wrrm = wx) -^cos- J * cosh ^ cos ? — i, (22) 

— oo 

where Ko((p) is the modified Bessel function of the second kind of argument <j>. In the CFT, 
where the correlation length tends to infinity, one is allowed to use the OPrH 

e(r)T(0) = ^ 2(Ak ~*~ Ar) O*(0), (23) 
k 

where Ok is a basis of fields, and C^j- are the dimensionless constants of the expansion. The 
most relevant operator appearing in (|23j) is the composite twist field Oo =: e7~ :, and it is defined 
implicitly as the twist operator which corresponds to the leading term. Its conformal weight 
and structure constant with e and T are known [8] 

A.. £T : = + A r and : = -. (24) 

In n 

We know that, due to the arbitrariness in the choice of the argument of Q : eT :, also derivatives 
of this field play a role in the OPE (f23j) . We are able to fix the weight and structure constants 
for these corrections, and to identify O a =: d 2a £T : with a £ N 

. 1 + 2a . , ~.a2a_T-. 1 

^~ +Ar and c " =wfv (25) 

where d 2a = {d 2 zl ) a . Notice that since e and T are both spinless operators only these particular 
derivatives can contribute to the expansion. From now on we will refer directly to the operators 
: d 2a £T :, denoting : eT : as the case a = 0. 

In the massive theory an analogue but different OPE exists [TT] which can be extracted 
peturbatively from (|23p . and can be used to express the LHS of (|22j) . Here and throughout this 
paper the assumption that there be the same basis of operators in the massless and massive 
theories is pushed forward, such that we will refer with the same symbols (with an abuse of 
notation) to the CFT and perturbed fields. This expansion is in a sense richer than the conformal 
one, and contains information about the scaling region sorrounding the conformal fixed point. 
This reflects on the OPE by giving a nonzero VEV to some of Ok, which are of nonlocal nature. 
In fact while in the CFT these values are identically zero due to scale invariance, in the massive 
case they are (Ok) ~ m 2Ak . The coefficients of this expansion contain information only about 
the short distance features of the theory, hence they are analytical with respect to the mass. 



*from now on the subscript labelling the manifold VEVs are evaluated on will be omitted. It is implicit that 
(T) is always evaluated on n disconnected copies of R 2 , while other VEVs are on the real plane. 

Hhis OPE is as usual to be taken in the "weak" sense; indeed it has a meaning only once plugged into a matrix 
element of the CFT on R 2 . 

^indeed the choice to take as argument is arbitrary, and one could have chosen any point in the interval 
[0, r]. The difference between these choices is represented by a Taylor expansion about 0, therefore the OPE (|23p 
should include all derivatives of the field : eT :. 
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Since in the Ising model the coupling constant is proportional to m, it is sensible to define an 
expansion of the form [11] 



C k eT (r) = cV(A fe -|-A r ) [l + C\mr + C\(mrf + ... 



By means of (|24p and (|25p we can express this expansion for the composite twist fields as 



ri:d 2a eT: 



'■■(r) 



l + 2a 1 

f 71 



n(2a)\ 



1 + C'f a£T '-mr + C'f" 6 ' 1 : (mr)' 2 + 



d 2a eT: 



for a G Nr 



(26) 



(27) 



Fixing the perturbation constants appearing in (J27]) is generally a very hard task already for the 
first order. However, due to the special nature of the OPE under consideration we are able to 
determine all C^ 2 " 6 ^ 1 in a systematic way. Therefore we expect the OPE in the massive theory 
to take the form 



e(r)T(0) = Y. C 'sT £T: (r) : d 2a eT : (0), 



(28) 



o=0 



although we see later that further corrections to (|28p . that are not predictable from CFT argu- 
ments, will also appear. 



3.1 Computation of (: d 2a eT :) and Cf 2 " eT: 

Let us now proceed by expanding the integral in the RHS of (|22[) in a small r region, and 
then compare the result with the massive OPE (j28|) . that is we compare terms with the same 
dimensions. In order to extract all the information needed it is convenient first to reexpress the 
integral in terms of the quantity t = mre? 



171 7T In 

—77 cos — mr » 
vr 2 2n v ' 



dt f 



Kn t + 



(mr)' 



V mr ) \ mr ) 



<-L.\ n -2C0S^ i±\ n +1 
V mr ) n \ mr I 



(29) 



One can notice that leading contributions to this integral are given for large t/(mr). This 
provides a natural parameter over which is possible to expand the fraction in (|29p as 



l-(-U~ 2 l-(-M~ 

\ mr J \ mr ) 



' t 
•. mr 



) " - 2C0S^ [-*-) » + 1 a 
/ n V mr ) 



^ Q (n) 



mr 



(30) 



where the coefficients a (n) are real numbers, evaluated in Appendix O Once ([30]) is plugged 
into ([29]) one gets 



— cos — Vf] a (n)( 
7r z 2n 



l+2a i 

mr ) n 



o=0 



l + 2a 



(it t — "]((, { + 



(mr)' 



(31) 



The dependence on mr of the Bessel function can be extracted by expanding it for small r as 
(mr) 2 



Kb t + 



K (t) 



-(mr) + 0(mr) 



(32) 



The resulting integrals can be thought of as a special case of a known integral of the Maijer G- 
function, as explained in Appendix [Bj Substituting (|32p into (|3ip and carrying out the integrals 
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gives 



< e(r )T(0)> = cos ^ f; 

a=0 



rt-l-2a\ 2 

J , ■ i 

2n y (mr) n 1 + 



+ r( ^7 } (mr)^+ n - ) r( ^7 } K)^% (33) 

ra+l + 2a 2 1+i: ^ V n + 1 + 2a J 2 2+2 ^ 

3 -p In— 1— 2a \ 2 

(n + l + 2a) 2 (3n + l + 2a) 2 2+2d ^ 

where we reported only the first contributions related to the composite twist field and its de- 
scendants, as we are mainly interested in those. After a bit of manipulation we can compare 
term by term this expansion with ([28j), and by matching the terms with the same perturbative 
order we are able to extract the following VEVs 

2a , _ cos (1+ 2 2 n a)7r (2a)!ra^ ( n - 1 - 2a\ 2 i+2a 



(■■ d 2a sT :> = r m^(T). (34) 

2 1+ « 7T 2 V 2n ) 

This is the main result of this section. In the same way we can fix the constants Cjd^eT: ^ Q 
every order. A major challenge when doing this is to be able, for terms proportional to the 
same power of r, to distinguish between contributions to expectation values and to structure 
constants. It turns out that this ambiguity can be resolved by requiring that the expectation 
values (|34p are continuous functions of n for each fixed value of a. This requirement is natural 
because of the special relation between T and the entanglement entropy, in which context it 
is necessary to analytically continue all physical quantities to n £ [1, oo). The first few non- 
vanishing coefficients are 

o , (l+2a)7r 



n+l + 2a (n + 1 + 2a) 2 tan ^~ 



2n 

T.d 2a eT: _ n2 n 3 (l + 2a + 2n) tan {1+ ^ )n 



C'f £l: =-, : ^77. : 7T^ + 



2(re + l + 2a)(3n + l + 2a) (1 + 2a + n) 2 (l + 2a + 3n) 2 tan £ 



2n 

(4n + 1 + 2a) 
) 2 (3n + 1 + 2a)(! 

(l+2a)7r „4/^ m , i , r, \ / 



n:d ia eT: _ 

6 6(ra + l + 2a) 2 (3n + l + 2a)(5n + l + 2a) 



sin^^ n 4 (4n + 1 + 2a) / _ 2n 



sin £ 2(n + 1 + 2a) 2 (3n + 1 + 2a) 2 (5n + 1 + 2a) V 1 + 3n + 2a, 

(35) 

We find that Cf^f'- = 0, for j 6 N . 

It is worth noticing that VEVs of these composite operators are singular whenever the argument 
of the Gamma-function is either zero or a negative integer. Analyzing ()34|) one can see that when 
n = l+2a the VEV of the 2a th derivative is divergent. In other words, such singularities can only 
occur for odd values of n. Therefore, it follows that the analysis above is only really consistent 
when restricting n to be even, in which case no singularities for special values of n and a arise. 
In the n odd case, the singularities that occur at various orders of the expansion actually cancel 
each other, generating a well defined short distance expansion. This stark contrast between 
the n even and n odd cases is a priori rather surprising. At present, we do not have a clear 
understanding of why it happens. It is a matter which we might came back to in future research. 
The leading term in (|28|) is well defined for all values of n > 1, for both n even and odd so that 
our identification of the expectation value of : eT :, already performed in [8], does still hold for 
general values of n. 
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3.2 Logarithmic corrections to the massive OPE 



Analyzing carefully the contributions reported in Appendix [Bj one can notice that further cor- 
rections involving terms of the form (mr) 2a and (mr) 2a log(mr), with a = 0, 1, . . . occur (note 
for example the logarithmic terms in (|83p and (|84j)). such that (|28[> is modified to 

oo 

e(r)T(0) = [ C ef a£T 'i r ) : ■■ (0) + mC^r T {r)d 2a T(0)] , (36) 

a=0 

where the new terms mentioned above are contained in the coefficients C^^(r). It is clear that 
they are features of the massive theory as they have no counterpart in the CFT. To understand 
their presence one has to consider once more (HJ). The term proportional to ip^ is responsible 
for all the composite fields : d 2a eT : in the OPE, whereas the term proportional to the identity 
generates contributions proportional to rriT and its derivatives md 2a T. It is interesting that 
this paves the way for the evaluation of (d 2a T), although in this case, due to the structure of the 
expansion, this would need a resummation of infinitely many terms. The presence of logarithmic 
terms is imputable to the freedom in choosing b in flip. Indeed in the expansions carried out in 
Appendix B we have general terms of the type 



oo 



(mr) 2a ( T («> A n ) + A («> 0. n ) lo g mr ) > ( 37 ) 

f3=0 

where T and A are two rational functions of a, (3 and n. When a correction of the kind of (|37p 
is plugged into ([55]) it gives for fixed a 

- cos —(mr) 2a V n 8 (n) (T(a, f3, n) + A(a, /3, n) log mr) . (38) 

13=0 

The term containing T contributes to (e)(d 2a T). The presence of the logarithmic term allows 
us to rewrite ([38]) as 

_ (DlH C o S l( mr ) 2a V ns(n) ( T(a, /?, n) + A(a, /?, n) log 5 + A(a, P, n) log ^f) , (39) 
7r z 2n \ o J 

13=0 

where 5 € M + \ {0}. This corresponds to a redefinition of (e), and one can notice that a 
logarithmic correction is the only functional form that allows this to happen. This is not the 
first time a logarithmic correction to two point functions of the Ising model has been observed 
(see for example the spin-spin correlation function in |15j). In both cases their presence is fully 
explained by the ambiguity in the definition of (e). 



4 Two particle form factors of composite operators 

This section is dedicated to the evaluation of form factors for composite twist fields of the 
type introduced in Section [3l These operators are formally defined as the regularised limit of 
an operator O approaching the twist field in the original CFT. The regularisation defines the 
meaning of the ordered product, which can be taken to be a point splitting 

: OT : (x) ~ lim 0(x + e)T(x). (40) 

One can start by considering : ipT ■ as a benchmark. Since T is even and ift is odd under the Z2 
symmetry of the Ising model, only odd particle form factors will be non- vanishing. Considering 
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then the matrix element (0| ip(x)T(0) \0) the one-particle form factor can be extracted by looking 
at the leading contribution as x approaches 0. It is useful to introduce the resolution of the 
identity for the n-copies system 



k=i m— /u/b 



d6i 



de, , 



e 1 >e 2 >->e k 



(27T)' 



L) • • • j Sk}m,...,fj, k n k ,...,ij,x\Ski ■ ■ • j &l\ 



(41) 



where /ij label the number of the copy and all quantum numbers of the i-th. particle. In this 
way one can write 



<O|#c)T(O)|0> 



n 

2n~ 



#<O|V(z)|0) <0|T(O)|fl). 



The matrix element (0| ip(x) \4>) can be easily extracted from © 

m(x)\4>) = J^£- ia *+, 

V 4-7T 

while (</>| T(0) \9) is linked to the two-particle form factor by crossing |12] 
(0| T(0) |0> = (0| T(0) |0, + in - ie + ) + (T>*(0 - <t>) 



(42) 



(43) 



(44) 



where the introduction of ie + has to be thought of in the distributions sense, and describes how 

to avoid the pole at m. 

Plugging ([MD and (@3) into (02J yields 



(O|^(s)T(O)|0) =(T)J^ef-^ + 

V 47T 



/' 



m /T\ 7F 

, — { / ) cos — 
2vr V 4vr w 2n 



sinh 



e 2 



In 



(45) 



sinh r ) sinh f 



2n y 



2« 



The leading term when x approaches is given by the integral part, that is divergent, but can 
be made convergent by shifting the domain of integration from the real axe K to M-Mvr/2. Such 
a change of variable does not affect the result of the integration, as the integrand has no poles 
in the region between the two axe^jj]. 
The integral in (|45p becomes then 



e '4 



sinh 



2ra 



sinh 



2n 



sinh 



-<j>-i3n/2 
2n 



(46) 



In this integral the main contributions come from large \(p\, so that, after splitting the integration 
path into positive and negative regions, one can find two series expansions for the fraction in 
l6|) . The leading term for small x is given by the positive ip part, that is 



-2e _i ?( 1+ ")e^ 
-2e~^ ( 1+ n) e ^ 
_2 e _i f ( 1+ n) e ^25( 1 "^)r 



(l-±)-x£, + _ 



e -mxcosh0 cosh 



^ I 1 - - 

n 



+ . 



(47) 



n- 1 
2n 



notice that it is the term ie + which allows to perform this shift. 
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The one particle form factor for : ipT : can be read out from (|45p and (|47p 




m2n 



\/4~7r 



e2n 



(48) 



The procedure used to get this result was introduced for descendants of twist fields in [13]. In 
that case the authors were dealing with the twist field associated to the global U(l) symmetry 
of the Dirac Lagrangian. In this and next sections extensive use of this procedure is made, 
demonstrating its consistency for a different kind of twist fields. 

4.1 Two particle form factor of : eT : 

One of the main focuses of this paper is the computation of form factors of composite operators 
defined as in (|40p with O = e. The matrix element (0| Q(x)T(0) \ 6162) can be written by means 
of gU as 



where j labels the copy number, and is repeated because e connects only particles on the same 
copy. 

The matrix element jj(<j>2, 0i| 7~(0) |0i, 62) is connected to the four-particle form factor of the 
twist field by the crossing relation (|44p . which used repeatedly gives 



jA<h<h\T(0) |0i0 2 ) = - <0| T(0) \9 1 ,9 2 , 0i + in - ie + ,4>2 +m- ie + ) lfljjtj + 
+6(6i - <h.)6i d (0| T(0) \6 2 , + ivr - ie + ) hj + 5(9 2 - 2 )Sij (0| T(0) \0 U 0i - in - ie + > 1>j - 
6{8i ~ <h)Sij (0| T(0) |0 2 , <Ai + »7r - ie + > 1J - <5(0 2 - fa)6 ltj (0| T(0) |^i, 2 - iir - ie + > l j + 
+6(0! - fa)5(9 2 - fa)Sij - 5(9 1 - fo)5(9 2 - fa)6 x>j . 



Even if (|50jl is quite cumbersome one can notice that many terms can be extracted from the 
others with the exchange 9\ < — > 92, which leaves the integration untouched. The leading 
contribution for small x is given by the first term (the one without deltas). This four-particle 
form factor considers two particles on the first copy, and two on copy j, but can be reduced 
to one where particles are considered on the same copy (say 1) with multiple applications of 
equations (fl"4"|) and (fl~5|) 



and this allows us to use the Pfaffian structure (|19p to re-express it. In the last step of (|5ip the 
notation (. . .)+ is introduced to indicate that any pole on the real axe of what is in the brackets 
is avoided with the ie + prescription. With the help of (|19p . (|20p and (|12p . one can rewrite the 
leading term of (|49|) as 




(50) 



<0|T(0) |0i,0 2 ,0i +Mr-ie + , 02 +in -ie + ) lXjij = 
Fp 11 " (9i,92, 0i + in - ie + ,4>2 + in - ie + ) = 

if' 1111 (01, 9 2 , 01 + (2j - 1)Z7T,0 2 + (2j - l)t7T) + 



(51) 




[F 2 ' |iX (0i,02K '"(01,02) - F 2 ' Vl,<h. + (2j - l)w)F 2 ' '"(02,02 + (2j - l)ivr)+ 
F 2 T|11 (0 2 ,0 1 + (2j - l)ivr)F 2 r|11 (0 1 ,0 2 + (2j - 



(52) 
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and proceed by evaluating the leading contribution for the three terms in (|52p . The first one is 



)S ^Fp 11 (9 l ,e 2 ) dfadfa , V 2 . { ) 2n J . , e-"K^), (53) 



sinh ( ^ J' 2 ) sinh 



in 2 " J sinh j S inh 

on which we can perform the change of variables t = (cfti — </> 2 )/2 and s = (<f>i + <j>2)/2, and carry 
out the s integration to get 

mcos-£- rm, , f sinh t sinh- , , , 

2tt 2 J smh (i + f ) smh (i - f ) 

We notice again that the leading contribution is given for large t so that we can use the parity 
of the integrand to reduce the region of integration to (0, oo), and expand the fraction in §J>4 
in a convergent way on this domain. The resulting integral is 

mcos 2^ P T|ii 

TT 2 



F 2 T|11 (0i,0 2 ) / dt e t ( 1 -k)K (mxe t ) + (55) 
J o 



and with the change of variable tt = mxe 1 we can extract the leading order for small x, that is 
_H^^ F T\ii i9id2) n dt u -hKo(u)Umx)^- 1 + .... (56) 
Solving the integral we finally obtain 



2 1 +£,r 2 V 2n 

The second and third terms in (|52p are more involved then the first, due to their explicit 
dependence on the parameter j. One can start by noticing that the value of the third term can 
be extracted from the value of the second one by changing sign, and performing the exchange 
9\ i — > # 2 . Hence we focus on the first term. As, in order to have a convergent integral, the 
integration axe has to be risen by iir/2, we need to study the pole structure of this term. The 
only kinematic poles which lie on the real axes arise from the cases j = 1, n and <fii = 9\ and 
<^2 — #2; but they are avoided with the ie + prescription. In general F 2 (9,4>+ (2j — l)«7r) has 
kinematic poles for <f> = 9 + 2(n — j + l)wr and <f> = 9 + 2{n — j)m. Considering that j runs from 
1 to n one can see that all poles group in even multiples of i7r, so that the first group above 
the real axe is in 9 + 2iir, and correspond to j = n, n — 1. Hence the needed shift can be safely 
performed. 

The integral we want to evaluate is then 

im(T) cos 2 n 



2 n 



87r 2 n 2 



/ #i#2 e-^i+^sinh 



sinh ( 01 ^ 2 ( W 2J 1>7r ) sinh , 



(58) 



sinh ( ei " 01 " 2 2 w (j " 1)CT ) sinh ( fr-y^ j sinh ( fe-*2-2(j-i)^ ginh ^ 2 -^-2ji« 

Although the fraction in (|58p is rather complicated and mixes integration variables with param- 
eters it can be dramatically simplified by means of the following identity 

sinh (a 1 — f3\ ± 7) sinh (a 2 — /3 2 ± 7) = 

1 (59) 
- [cosh (ai + a 2 - (/?i + /3 2 ) ± 27) - cosh (ai - a 2 - - /3 2 )] , 
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leading to 

-.2 7T n 



cQsh ^ 1 +O 2 -(0 1 +^)-2(2j-l)i^ _ cogh ^^_^-M 

cosh (M^ ^H^ - cosh ^i-^-^i^ T) 

1 

cosh f gL -^-to+feMM _ cosh ( SgE ^T 



(60) 



Now performing the same change of variable as in (|54j) gives 



im(T) cos 



2 jr_ « 

2n 



^ J dtds e -2ir^sinh S co S ht sinht 



(61) 



2vr 2 n 2 

3=1 

cosh (^-^^)- cosh 
cosh _ £±2^) _ CO sh (fi^-i 

1 

cosh _ ^) _ cosh li 

and to make it convergent the shift s — > s — in/ 2 has been performed. 

As x approaches the main contribution is given by large t, and s peaked around 0. It is natural 
then to expand the fraction in (|6ip in powers of t. As before there is no such expansion on the 
whole real axe, but splitting it into the positive and negative parts, allows us to consider two 
different series which converge respectively on the two regions. We can start considering t < 0, 
such that the expansion of (|6ip yields 



oo 



2e~" f° dt e-'( 1 -^)Ao( 

J — oo 



(62) 



mxe ) + .... 



Following now the same procedure which was used to obtain ()57p out of (|55p . we arrive at the 
result 

im(T) cos 2 /ra-l\ 2 ei-o 2 ., J -. „ s 

2n T — — e — ^H"" 1 + .... (63) 



2- +1 vr 2 n V 2re 

The positive part of the integral can be performed with the same logic and gives 

im(T) cos 2 (n - 1\ 2 »i-9 2 . , i , . 

!n T — — e^r(mi)" 1 + ..., (64) 



2^+Vn V 2n 
so that the final result of the second part is 

im (T) cos 2 § p /n-l\ 2 , ^ 9i-9 2 / i.. 

As mentioned before, the result of of third part of the integral in (|52|) can be obtained from (|65|) 
by switching 6>i « — >• 02 with a minus sign in front, which doubles the result. Putting ([57]) . ([63]) 
and ()65p together, we finally obtain the two particle form factor for the field : eT :, that is 

2 



cos £ / n - 1 



2 1 +t 7 r 2 V 2n 



ifl^AH 1^ (T ) sinh^ 



n 2re 



(66) 
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This result was the aim of this section. To check its validity one can employ (|16|) . Indeed, since 
: eT : is still a twist field, the same type of residue equations as for T must be satisfied. Then 
one can easily check using ([34"]) and ([66]) that 

lim(0 - 6)F' 2 £TAll {9 + wr, 9) = i(: eT :), (67) 
e->e 

which confirms the compatibility of the two results of these sections. 

Notice that this result satisfies all form factors equations for the twist field, and has a structure 
of the type 

if T:|11 (0i,0 2 ) = a [Ql(0i,e 2 ) + Pk(6 1 ,6 2 )] F min (9 1 ,e 2 ), (68) 

where a and (3 are two dimensional constants, F m i n is the minimal form factor of the theory!**!, 
and k is a kernel solution of the form factor equations (see e.g. [5] for a discussion). Even if we 
could have understood that ([66]) should have had the form ([68]) to fulfill all twist properties it 
has to, we would have never been able to fix a and f3 without the methods used in this section, 
and Section [3] 



5 Higher particle form factors of composite operators 

In this section we deal with the computation of higher particle form factors for : ipT : and : eT : . 
They can be extracted with the same methods for both operators, employing higher particle 
form factors of the twist field (|19j) . Indeed, focusing on : ipT :, when looking for the leading 
term for the 2k — 1 particle form factor, one has to deal with 



<0| ^(x)T(O) 



n 



>2k-\i 



— I #<0|#z)|tf>><0|T(0) 



hk-x, 4> + wr - . 



One is then able to isolate term by term the higher particle part, and reduce it to the same 
evaluation carried out in (|42 |) -(|48 |) . getting finally 



F, 



:i/>7~:|ll...l 



2k-l 



■>pT:> 



(70) 



where K-^q-- is the 2k x 2k matrix defined as 



K. 



ipT: 







F 



\- F -f T --\\9 2k . 1 ) —F 2 T][11 (d\,Q 2 k-i)/{T) 



Fi i,T --\\e 2k ^) 



? T\ll 



(0 



1) V2k-1 







)/{T) 



(71) 



/ 



Higher particles form factors for : eT : can be evaluated with the same logic, although they 
show a more complicated pattern, and can not be reduced to a Pfaffian form. This is due to the 
presence of a kernel part in the two particle form factor. The 2k + 2 particle form factor is 



7 : £ T:|ll...l, 

2k+2 \ 



7 1,V2, 



>2k+2\ 



E 

Kj 



4 Wm) 



<r> 



2k ± 2k 



h, - ■ ■ ,02k+2)ijT 2 



:eT:\ll 



(72) 



where cr(i,j) is the permutation that brings 0i and 8j to the right of all other rapidities, while 
with F(. . . )ij we mean a form factor of all rapidities but those two. 



**for the Ising model -Fmin (Si, ^2) = — isinh^. 
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6 Conclusions 



In this paper we have investigated the short-distance behaviour of the correlation function 
(e(r)T(O)) for the two dimensional Ising model in the vicinity of the critical point beyond the 
leading contribution. This led us to the identification of the vacuum expectation values of a new 
class of twist fields, including 



and its derivatives. Furthermore we managed to compute massive corrections to the structure 
constants up to (mr) 6 for all these fields. The very fact that we are able to compute VEVs of 
such a large set of local operators is remarkable, as the computation of VEVs for general theories 
and fields is known to be a very hard task and there is no general procedure to tackle it. In 
addition to the contributions that would be naturally expected from the CFT theory we have 
found logarithmic corrections, a phenomenon that has already been observed for other correla- 
tion functions in the off-critical Ising model (see e.g. [I5|). These logarithmic terms are due to 
the arbitrariness in the definition of (e). In addition we have computed all higher-particle form 
factors of : eT : and : ipT ■■ By exploiting (|40p we have been able to fully determine the nor- 
malization of all form factors and to provide new solutions to the form factor equations for twist 
fields. A byproduct of our investigation is the fact that all the expectation values (: d 2a eT :) 
are negative. Although we did not put much emphasis on this feature in the present paper, this 
provides further evidence for the negativity of the connected correlator (e(r)T(O)), claimed in [8J. 

There are a number of open problems related to the present work which we would like to 
address in future. First, the short distance expansion of e(r)T(0) that we considered here is 
only well defined for even n. It would be interesting to have a deeper understanding of what 
happens when we consider n odd, and why these two cases are distinct. Second, this paper 
introduces an OPE which involves the twist field in a replica theory. This type of OPEs have 
never been studied, neither in the massive nor in the critical theories. A better understanding of 
the operator content and OPEs in replica theories is desirable. Finally, it would be interesting 
to understand if these composite twist fields are in any way related to the entanglement entropy 
of particular states in the Ising model. 

Acknowledgements The author wishes to thank Olalla Castro-Alvaredo for her guidance 
and for proofreading this manuscript (without her help the completion of this project would not 
have been possible), and Benjamin Doyon for precious explanations and discussions. 



A Coefficients Q a (n) 

This section is devoted to the computation of coefficients in ()30|) . As explained in Section [3] the 
aim is to expand the fraction on the LHS of that equation for large t/mr. First of all, let us 
introduce the more convenient variable u = (t/mr)~ 2 / n . The denominator can be then treated, 
as long as t > mr as the generating function of the Chebishev polynomials of the second kind, 
that is 



for — 1 < x < 1, and |u| < 1. The first condition is satisfied for every n > 2 as x = cos(7r/n), 
while the second is satisfied in the whole integration domain of (|29p except for the lower limit 
t = mr. This divergence is "cured" by integrating over the domain [mr + e, oo), where e is a 



: eT : (x) ~ lim e(x + 5)T(x) 



(73) 




(74) 
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small parameter. Once this expansion is plugged into (|29p the sum and the integration can be 
safely exchanged. After the integration is performed one has then to be sure that the result does 
not depend on e, and finally set it to zero. We have performed these steps showing that indeed 
()74p in this case can be taken as valid also at the point t = mr. The details are technical and 
cumbersome, and they will not be reported here. From now on, and throughout the calculation 
in Section [3] we will take (|74|) as series representation an the whole integration path. 
The polynomials U a (x) in our case are formally defined as follows 

■ (1+«)tt 

/ 7T \ sin j — 

U a (coB-)= . I , 75 
V n J sin - 

n 



The LHS in (1301) can then be expanded as shown in the RHS with 



£l a (n) 



cos J — — — . r 

— — 4r — if a < n 

(l + 2a)n . (l + 2aW • \ IU ) 

COS J — -= — Sin i — 7= — . „ 

¥ 1 — if a > n 



B Definite integrals of Bessel functions and powers 

In this appendix we present a solution to integrals of the kind 

/•oo 

/ dt t-»K v {t), 

Jmr 



(77) 



where both \i and mr are positive real numbers. In Section [3l in particular, an expansion for 
small values of mr was needed, so that this will be the aim of this appendix. First of all let us 
introduce the function 



si m,n 
P.9 



ai, 



, Up 



h b 



) "q 



1 

2tH 



n7=i r (^-^)n- = ir(i-a J + S ) 



rr<? 
L l lj=m+ 



i r(i-6 i + s )n p =re+1 r(a i - s ) 



t s ds. 



(78) 



This is a representation of the Meijer G-function, in the formalism adopted by |T3], and the 
details and properties about this function will not be reported here. A useful identity is 



Ku{t) = \ G 2 ; 2 ° 



v —v 
2'1T 



(79) 



which holds for | arg(t)| < ir/2, and the empty sets of gamma functions' poles are omitted. In 
the light of (fT9|) the integral in ([77]) can be rewritten as follows 



(mr) 



1 — poo 



dt t — 5- Gq 5 1 —t 



V V 

2'~2 



(80) 



This is a special case of a known integral of the G-function, which in the most general form is 



/oo 
dt t-P{t-Vf- x G%* 



at 



h,...,b q 



T(a)Gl 



m+l,n 



O: 



Ol, . . . , a p , p 



P ~ er, 61 



,...,b q 



(81) 



which holds for real | arg(t)| < (m + n — p/2 — q/2)ir, p + q < 2{m + n), K(<t) > and 
!R(p — a — dj) > —1 Vj G [l,ra]. These conditions are all clearly satisfied by (|80p . so that the 
result is 



dt t-^Kyit) 



(mr) 1 M „3o/ m2 



2 

fi—l u_ _u_ 
2 ' 2 ' 2 



(82) 
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Now this result has to be restricted to the cases (|32[) to be useful for the OPE that was considered 
in Section [3l The calculations are tedious and the results cumbersome, hence only the first few 
terms of the expansion of the first two contributions are given 

°° , 2a + 1 t » 1 l + 2q _ ( TL — 1 — 2tt \ 2 

A r— K f =2- 1 -— r — + 

V 2n ) 

n\(n — 1 - 2a) (7 — log 2) — n] . , i 1+2^ n , ,-, i+2 Q 

+ ^ (n _i_ 2a)2 W - + (re _ 1 _ 2a) (^) 1 - I0g(mr)+ 

n[(3n — 1 — 2a) (7 — 1 — log 2) — n] , .0 i+2a n . i+2a 

+ 4(3n- l-2a)3 (mr) " + 4(3n - 1 - 2a) ^ ^ bg(mr)+ 

. c l + 2a 

+ 0{mrf n 

(83) 



r a r 1 -^^) = (mrrl " 1 ''" re + 2- 2 -^r f n " 1 " 2a l r f- 1 + 2a + n l + 

7 mr iW l + 2a + n V 2n J \ 2n J 

n[(n — 1 — 2a)(l + 21og2 — 27) — 2nl 1 1+2^ n . , ■■ i+2a 

4(n - 1 - 2a)^ (mr) " " ~ 2(n - 1 - 2a) (mr) " " l0g(mr)+ 

n\(3n — 1 — 2a) (5 + 4 log 2 — 47) + 4nl . .0 i+2a n , l±2a , 

+ 64(3n-l-2a)' V) 3 "— + ^ _ i _ 2q) log(mr)+ 

r l + 2a 

+ 0{mrf 

(84) 

where 7 = 0.577216 is the Euler-Mascheroni constant. 



References 

[1] V.G. Knizhnik, Analytic fields on Riemann Surfaces II, Comm Math Phys 112, 567-590 
(1987). 

[2] P. Calabrese and J. L. Cardy, Entanglement entropy and quantum field theory, J. Stat. 
Mech. 0406:P06002 (2004). 

[3] O. A. Castro-Alvaredo and B. Doyon, Bi-partite entanglement entropy in massive 1+1- 
dimensional quantum field theories, J. Phys. A 42, 504006 (2009). 

[4] B. Doyon, Bi-partite entanglement entropy in massive two-dimensional quantum field the- 
ory, Phys. Rev. Lett. 102:031602 (2009). 

[5] O. A. Castro-Alvaredo and E. Levi, Higher particle form factors of branch point twist fields 
in integrable quantum field theories, J. Phys. A 44, 255401 (2011). 

[6] G. Delfino, P. Simonetti and J. L. Cardy, Asymptotic factorisation of form factors in 
two-dimensional quantum, field theory, Phys. Lett. B387, 327-333 (1996). 

[7] A. B. Zamolodchikov, Irreversibility of the flux of the renormalization group in a 2-D, field 
theory, JETP Lett. 43, 730-732 (1986). 

[8] O. A. Castro-Alvaredo, B. Doyon and E. Levi, Arguments towards a c-theorem from 
branch-point twist fields, J. Phys. A 44, 492003 (2011). 



15 



[9] A. B. Zamolodchikov, Integrable field theory from conformal field theory, Adv. Stud. Pure 
Math. 19, 641-674 (1989). 

[10] J. L. Cardy, O. A. Castro-Alvaredo and B. Doyon, Form factors of branch-point twist 
fields in quantum integrable models and entanglement entropy, J. Stat. Phys. 130, 129-168 
(2008). 

[11] Al. B. Zamolodchikov, Two point correlation function in scaling Lee- Yang model, Nucl. 
Phys., B348, 619-641 (1991). 

[12] F. A. Smirnov, Form-factors in completely integrable models of quantum field theory, Adv. 
Ser. Math. Phys. 14 (1992). 

[13] B. Doyon and J. Silk, Correlation functions of twist fields from Ward identities in the 
massive Dirac theory, J. Phys. A 44, 295402 (2011). 

[14] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series, and Products, Academic Press, 
New York and London, 1965. 

[15] B. Mikhak and A. M. Zarkesh, Manifestly Finite Perturbation Theory for the Short- 
Distance Expansion of Correlation Functions in the Two Dimensional Ising Model, 
Nucl.Phys. B430, 656-682 (1994). 



16 



